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^D . Abstract 

(N ■ 

We study the geodesic equation in the space-time of an Abelian-Higgs string and discuss the motion 
of massless and massive test particles. The geodesies can be classified according to the particles energy, 
angular momentum and linear momentum along the string axis. We observe that bound orbits of massive 
particles are only possible if the Higgs boson mass is smaller than the gauge boson mass, while massless 
particles always move on escape orbits. Moreover, neither massive nor massless particles can ever reach the 
string axis for non-vanishing angular momentum. We also discuss the dependence of light deflection by a 
cosmic string as well as the perihelion shift of bound orbits of massive particles on the ratio between Higgs 
and gauge boson mass and the ratio between symmetry breaking scale and Planck mass, respectively. 



o\ 



^ 1 Introduction 

en " 

;q Cosmic strings have gained a lot of renewed interest over the past years due to their possible connection to string 
00 . theory pQ . These are topological defects [5] that could have formed in one of the numerous phase transitions in 
^^ ' the early universe due to the Kibble mechanism. Inflationary models resulting from string theory (e.g. brane 
[~^. . inflation) predict the formation of cosmic string networks at the end of inflation [3] . It would be very interesting 
< T^ ' to observe these objects in the universe. In recent years the detection of these objects has focused on the Cosmic 
^D , Microwave background (CMB) data |3], though other observational effects such as gravitational lensing have 
T 7". ' also been discussed [5]. Here we discuss the possibility that cosmic strings might be detected due to the way 
^ . that test particles move in their space-time. 

Until now the geodesic motion of test particles in space-times containing cosmic strings has been mostly 
studied in the limit of vanishing width of the cosmic string. Since the space-time of an infinitely thin cosmic 
Cd ' string is locally flat [2] geodesies are just straight lines. In [6] a general cosmic space-time has been studied and 
it has been shown that the geodesies of massless particles must move to infinity in both direction. Geodesies in 
cosmic string space-times have also been used to explain the motion of particles in elastic solids [7j . Furthermore, 
black hole space-times containing infinitely thin cosmic strings have been investigated both for static black holes 
8, 9, [10l [Tl] as well as for rotating black holes (9j H21 [13l HH [15] . In this case, analytic solutions to the geodesic 
equation in terms of elliptic functions are possible. 

In this paper we are aiming at understanding particle motion in the space-time of a finite width cosmic string. 
The underlying field theoretical model is the U(l) Abelian-Higgs model which has string-like solutions [T5] . 
The gravitational properties of Abelian-Higgs strings have also been studied in detail by minimally coupling 
the Abelian-Higgs model to gravity. Far away from the core of the string, the space-time has a deficit angle, i.e. 
corresponds to Minkowski space-time minus a wedge |17j . The deficit angle is in linear order proportional to 
the energy per unit length of the string. If the vacuum expectation value (vev) of the Higgs field is sufficiently 
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large (corresponding to very heavy strings that have formed at energies much bigger than the GUT scale), 
the deficit angle becomes larger than 2tt. These solutions are the so-called "supermassive strings" studied in 
[18| and possess a singularity at a maximal value of the radial coordinate, at which the angular part of the 
metric vanishes. Interestingly, it was realized only a few years ago [19l [20] that both the globally regular 
string solution as well as the supermassive solution have "shadow" solutions that exist for the same parameter 
values. The string-like solutions with deficit angle < 2tt have a shadow solution in the form of so-called Melvin 
solutions, which have a different asymptotic behaviour of the metric and have higher energies than their string 
counterparts (and are thus very likely cosmologically not relevant). The supermassive string solutions on the 
other hand have shadow solutions of Kasner-type. Kasner solutions possess also a singularity at some finite 
radial distance, the difference is that for Kasner solutions, the it-component of the metric vanishes at this 
maximal radial distance, while the angular part of the metric diverges there. Since both supermassive as well as 
Kasner solutions contain space-time singularities they are surely of limited interest for cosmological applications. 
In this paper we thus concentrate on the motion of test particles in the space-time of a gravitating cosmic string 
that possesses a deficit angle < 2ir. Since the field theoretical solutions and in particular the metric functions 
can only be determined numerically analytic results are not possible for the geodesic equation. This is similar 
to the geodesic motion in the space-time of a gravitating magnetic monopole 1211 . where also only numerical 
results are possible. 

Our paper is organised as follows: in Section 2, we discuss the field theoretical model that possesses string- 
like solutions and we also give the geodesic equation. In Section 3 we discuss our numerical results, in particular 
we give examples of orbits and demonstrate how the ratio between the symmetry breaking scale as well as the 
ratio between the Higgs and gauge boson mass influence our results. We conclude in Section 4. 

2 The model 

In the following we will first present the space-time of an Abelian-Higgs string and then give the geodesic 
equation that describes the motion of massless and massive test particles in the space-time of an Abelian-Higgs 
string. 

2.1 The space— time of an Abelian-Higgs string 

The Abelian-Higgs model is a field theoretical model with cosmic string solutions [16]. Here, we couple this 
model minimally to gravity. The action then reads 

where R is the Rlcci scalar and G denotes Newton's constant. The matter Lagrangian C m is given by 

C m = D^(D^)* - \f^F^ - A (#* - iff (2) 

with the covariant derivative D^(j) = V M — ieA^cj) and the field strength tensor F^ v — d^A v — d^A^ of the U(l) 
gauge potential A^ with coupling constant e. The field is a complex scalar field (Higgs field). 

The most general, cylindrically symmetric line element invariant under boosts along the z— direction is 

ds 2 = N 2 (p)dt 2 - dp 2 - L 2 {p)d<p 2 - N 2 {p)dz 2 . (3) 

For the matter and gauge fields, we have [16] : 

4>(p, V ) = vh(p)e inv , A lt dx* = ±(n-P(p))dtp, (4) 

where n is an integer indexing the vorticity of the Higgs field around the z— axis. As is apparent from the 
Ansatz the cosmic string possesses a magnetic field along the z-axis with [TS] 

Bzip) = - eL(p) dp ■ (5) 
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There are three mass scales in our model: the Higgs mass M H = vAry, the gauge boson mass M w — \/2erj 
as well as the Planck mass Mpi = G -1 / 2 . The Bogomolny limit A = 2e 2 [22] corresponds to equal gauge and 
Higgs boson mass. The cosmic string has finite width. In fact there is a scalar core with width /?h 
a magnetic flux tube core with width pw ~ A^w ■ ^ n ^ ne Bogomolny limit pn = Pw- 

We can then do the following rescaling 

P T L 

er\ er\ 

such that the total Lagrangian only depends on the following dimensionless coupling constants 

A 



7 = 8ttG^ , /3 = 



(G) 



(7) 



Note that 7 is proportional to the square of the ratio of the symmetry breaking scale r\ and the Planck mass 
A/pi, while f3 is proportional to the square of the ratio between the Higgs boson mass Afn and the gauge boson 
mass Afw 

Varying the action with respect to the matter fields and metric functions, we obtain a system of four 
non-linear differential equations. The Euler-Lagrange equations for the matter field functions read: 
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while the Einstein equations are 
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Here and in the following the prime denotes the derivative with respect to p. 

The set of differential equations can only be solved numerically subject to a set of boundary conditions. The 
requirement of regularity at the origin leads to the following boundary conditions 



h(0) = , P(0) = n , N(0) = 1, iV'(O) = 0, L(0) = , L'(0) = 1 , 
while the finiteness of the energy per unit length requires 

h(oo) = 1 , P(oo) =0 . 



(12) 
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In this paper we are interested in test particle motion in the space-time of a cosmic string with deficit angle 
A < 27r, i.e. we do not consider supermassive strings or space-times of Melvin- or Kasner-type in this paper. 
The metric functions then have the following behaviour at infinity 



iV(p>l) = ci , L(p> 1) =c 2 p + c 3 



(14) 



where c\ , ci and C3 are constants that depend on /3 and 7 [20] . We have c\ = 1 independent on 7 and C2 = 1 — 27 
in the Bogomolny limit /3 — 2. c\ and c 2 both decrease for fixed 7 and increasing /3 with ci > 1, C2 > 1 — 27 
for /? < 2, and C\ < 1, c% < 1 — 27 for /3 > 2. Moreover, c\ increases with 7 for fixed j3 < 2, while c\ decreases 
with 7 for fixed /3 > 2. 



In the Bogomolny limit j3 — 2 which corresponds to Mh = M-yv we have that P' = L(h 2 — 1) [22]. Inserting 
this into the Einstein equation (ITU1) and using the boundary conditions it is easy to see that in this case N(p) = 1, 
i.e. the space-time is only determined by L(p). 

We define as inertial mass per unit length of the solution 

M = / V^93T$dpd<p (15) 

where <?3 is the determinant of the 2 + 1-dimensional space-time given by (t, p, (p). This then reads: 

oo 

p = 2kJnL ([h'f + ££ + ^! + I (h 2 - I) 2 ) dp (16) 

o 

In flat space-time 7 = and in the Bogomolny limit /3 — 2 the energy per unit length is directly proportional 
to the winding number n and is given hy p — 2irn |22j , while the general expression for p in flat space-time is 

p = 2nnf(/3) (17) 

where / is a slowly varying function of /3 which is / < 1 for < 2, / = 1 for ft = 2 and / > 1 for f3 > 2. For 
7 > and fixed /3 the mass per unit length of the solution decreases with increasing 7. 

In linear order the deficit angle A of the space-time is proportional to the product of the gravitational 
coupling and the energy per unit length of the string 7/i. Taking the non-linear effects into account A is related 
to c 2 (see (PJ|)) by 

A = 2rr(l-ca). (18) 

Following the discussion above, the deficit angle increases when increasing 7 and A = An^y for /3 = 2, while 
A < 47T7 for /3 < 2 and A > Att^ for j3 > 2. Since we are only considering space-times with A < 27r in this 
paper there are restrictions on the parameters 7 and j3. In [2D] the value of 7 at which A = 2n has already been 
investigated for some particular values of f3. We come back to this in the numerical results section. 

2.2 The geodesic equation 

We consider the geodesic equation 

where T^. denotes the Christoffel symbol given by 

K° = \^ ( 9 P9™ + d °9 P v - d v9pa ) (20) 

and t is an affine parameter such that for time-like geodesies dr 2 = g^dx^dx" corresponds to proper time. 
The geodesic Lagrangian C g for a point particle in the space-time ([3]) reads 
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where e = for massless particles and e = 1 for massive particles, respectively. 

The constants of motion are the energy E, the angular momentum of the particle that is aligned with the 
axis of the string (here the z-axis) L z and the momentum p z in z-direction 

„ , ,0 dt n dtp ndz , 

E-=N 2 — , L z :=L 2 -f , Pz :=N 2 —. 22 

dr dr dr 



Using these constants of motion and the rescaling ^ as well as letting £ — > £/{e 2 rj 2 ), p z — > p z /(er]), L z — > 
L z /(e 2 r] 2 ) we can rewrite (f2"Tj) as follows 

~(¥) 2 =£-V eS (p), (23) 



2 \dr r 
where £ = (E 2 — e)/2 and V c g (p) is the effective potential 
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As is obvious from (j2"3")l test particle motion is only possible for £ — V c s (p) > 0. The p-component of the geodesic 
equation ^§ = -^j* then reads 

For our numerical calculations, we have rewritten the components of the geodesic equation in the following way 

d V = ± -4^-r. —^ ( 26 ) 
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dz = ± -jjz , (27) 

dt ± — EdP ~ TT/2 ■ ( 28 ) 
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For positive angular momentum L z the positive and the negative sign denote outward and inward motion, 
respectively. The motion becomes non-planar for p z =/= 0. 

3 Numerical results 

The equations of motion (|8|)- (111[) can only be solved numerically. We have done this using the ordinary dif- 
ferential equation solver COLSYS [23 . The relative errors of the numerical integration are on the order of 
10~ 10 — 10~ 13 . The numerical data for the metric functions N and L then has to be interpolated. This has 
been done with a piecewise cubic Hermite interpolating polynomial within MATLAB. The interpolated data 
can then be used to find t(p), <p{p) and z(p) via (f2l>l) - (f2"5)) . This latter integration has been done by using the 
recursive adaptive Simpson quadrature within MATLAB requiring an absolute error tolerance of 10 -8 . 

As mentioned above we are only interested in cosmic string space-times with A < 2-k in this paper. Since 
A depends on both 7 and /3 there are restrictions on these parameters. These restrictions have already been 
investigated in [20] for some particular values of f3. Here, we show the domain of existence of cosmic string 
solutions with A < 2tt in FigQ] where we give the value of 7 max , i-e. the value of 7 at which the deficit angle 
becomes equal to 2n in dependence on j3. 

In the following we will distinguish between bound orbits and escape orbits. Note that when we talk about 
bound and escape orbits we are referring to the motion in the x-y-plane. The particles can, of course, move 
along the full z-axis from —00 to +00 for p z ^ 0. 

Bound orbits are orbits on which test particles move from a minimal value of p, p mm > to a maximal value 
of p, Pmax < 00 and back again. These orbits have hence two turning points with (dp/dr) 2 = 0. On escape 
orbits, on the other hand, particles come from p = 00, reach a minimal value of p, p mm > and move back to 
p = 00, which means that escape orbits have only one turning point with (dp/dr) 2 = 0. Looking at (|23j) it is 
obvious that turning points are located at those p at which £ — V c ff(p) = 0. 

For all our numerical calculations we have chosen n = 1. 
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Figure 1: We show the value of 7 max , i.e. the value of 7 at which the deficit angle A is equal to 2-k in dependence 
on f3. Below this line cosmic string solutions with A < 27T exist, while above this line the cosmic string solutions 
are supermassive with A > 2tt and hence possess a space-time singularity. In this paper we are only interested 
in solutions with A < 2ir. We also show 7 cr , i.e. the value of 7 such that for 7 > 7 cr the Gaussian curvature K 
of the 2-manifold with coordinates (p, (p) can have negative values for some p. 



3.1 The effective potential 

The first observation is that since N(p — > 0) — > 1 and L(p — > 0) — > p there is always an infinite potential barrier 
at p — > for L z 7^ which means that the particles can never reach the z-axis, i.e. p — 0. For L z — this 
potential barrier disappears. 

Moreover, in the Bogomolny limit (3 = 2 with N(p) = 1 the effective potential is monotonically decreasing 
from V c ff(p = 0) = +00 to V c ff(p = 00) = and £ — V e g(p) = can only be fulfilled for one p. Hence in the 
Bogomolny limit bound orbits do not exist and test particles can only move on escape orbits. This is different 
for [3 y^ 2 where N(p) is non-constant. In fact, all our numerical results indicate that bound orbits exist only 
for (3 < 2 and for massive particles (e — 1). An indication that this is correct can be seen when looking at ([23| . 
Bound orbits are possible whenever V e e possesses at least one extremum, i.e. there is one finite p for which 
dVefi/dp = 0. Using the explicit form of the effective potential (|24|) this gives 
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Assuming N and L to be positive, the right-hand side of (J29I) is always positive. In addition L' is always positive 
which can be seen from (|I1[) and the fact that L'(0) = and L'(oo) > 0. Now from (|1II) it is easy to see that 
V = for some finite < p < 00 is excluded. Hence, for extrema to exist we need to require that N' > 0. 
In fact for j3 < 2 we find that N' > 0, while N' < for /3 > 2 (see also the results in PH [2"0]), Of course, 
the requirement that N' > does not guarantee that bound orbits exist (this still depends on the choice of £, 
L z and p z ), but for N' < 0, i.e. j3 > 2 we can exclude the possibility of bound orbits. Our numerical results 
demonstrating the change of the effective potential with 7 and (3 are shown in Fig(2] For fixed 7 it is obvious 
that local extrema exist for (3 < 2, while there are no extrema for (3 > 2. Moreover, the potential changes only 
little for changing (3 > 2. 

For fixed j3 < 2 (here j3 = 1.25) local extrema do exist and the value of these extrema is enhanced with 
increasing gravitational coupling. For fixed (3 > 2 (here (3 = 2.5) no extrema exist. While the effective potential 
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Figure 2: The effective potential for the motion of a massive test particle (e — 1) in the space-time of an 
Abelian-Higgs string is shown for fixed 7 = 0.3, E — 1.009, p z = 0, L z = 0.03 and different values of (3 (a), 
fixed /3 = 1.25 < 2, E = 1.029, p z =0, L z = 0.1 and different values of 7 (b) and fixed j3 = 2.5 > 2, JS = 1.004, 
p 2 = 0, L 2 = 0.02 and different values of 7 (c). The black horizontal line represents the value off. 



possesses local extrema for /? < 2 bound orbits exist only for particular choices of S , L z and p z . As mentioned 
above, our results indicate that bound orbits exist only for massive particles (e — I). In order to demonstrate 
this, we have chosen j3 = 0.5, 7 = 0.4. We show the effective potential for a massive and a massless test particle 
in FigEJand FigHl respectively. It can be seen from Fig|3]that the potential possesses local extrema and that 
£ intersects the potential twice for sufficiently low values of E if e = 1. This is different for the massless case. 
As can been seen from FigfJ] the local extrema of the effective potential disappear with decreasing E such that 
when £ is comparable to the asymptotic value of the effective potential, there are no local extrema at all. Hence, 
all our numerical results indicate that bound orbits for massless particles are excluded. 

In the following we want to compare this result with the analytical result of [6], which states that for a 
general cosmic string space-time with topology R 2 x E massless test particles must move on geodesies that 
escape to infinity in both directions, i.e. closed geodesies are not possible. The assumption made in [6] is that E 
must have positive Gaussian curvature. To show that E has positive Gaussian curvature in our case, we rewrite 
the metric ([3]) for massless particles (ds 2 = 0) moving in a plane parallel to the x-y-plane as follows 
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dt = —dp + -r^dif =g lj dx l dx :l , i = l,2 



(30) 



where <jij is the so-called optical metric [23] of which the spatial projection of geodesies of massless particles, 
i.e. light rays are geodesies, c/ij is the metric of the above mentioned 2-manifold E and has Gaussian curvature 
K given by 
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where in the second last equality we have used the Einstein equations (fTO)) . ([TT]) with Tq the energy density and 
T£ the pressure in ^-direction. Note that 



made in that paper, i.e. 
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([3Tj) reduces to the relation found in [6] if we use the assumptions 

El. 



In the BPS limit we know that N = 1 and the Gaussian curvature is obviously positive, away from the BPS 
limit one has to use the numerical solution and compute the curvature. We find that for most values of /3 and 
7 the Gaussian curvature is indeed positive and our result is in agreement with that of [BJ. However, if the 
ratio between Higgs and gauge boson mass is sufficiently large, we find that K can become negative close to the 
string axis. Our results are shown in Fig|TJ where we present 7 cr in dependence on j3. For 7 > 7 cr the Gaussian 
curvature K can become negative for values of p close to zero, while K is strictly positive for 7 < 7 cr . Though 
the theorem of [6] is not applicable here, we nevertheless find that bound orbits do not exist. 
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Figure 3: The effective potential for the motion of a massive test particle (e = 
Abelian-Higgs string with /3 = 0.50 and 7 = 0.40 for different values of E. Here p z 
red horizontal line represents the value of £. 



1) in the space-time of an 
= 0.100 and L z = 0.18. The 



In Fig|5] we show how the effective potential changes with changing p z and L z . Increasing p z shifts the 
effective potential simply to higher values which is apparent when looking at (|24|). Moreover, the lower the 
value of L z the lower the value of the minimum of the potential. Moreover, the potential barrier at p = 
disappears for L z = 0. 

We have then investigated how the existence of escape and bound orbits depends on the choice of E, L z and 
p z . Following the known discussion of orbits in the Schwarzschild case [25., we define v := L 2 Z and p :— E 2 . In 
FigEJa) we show the p-v-p\a.ne of a massive test particle for /? = 0.5, 7 = 0.4 and p z = and indicate for which 
values of p and v bound orbits exist. This is only possible in regions M2 and M 3. The corresponding effective 
potential V c g for regions M1-M4 is shown in FigJTJ In region Ml the effective potential V c g is always larger 
than £, so no particle motion is possible. In region M2 the potential has a local minimum and £ intersects V e g 
twice. The intersection points correspond to the minimal and maximal radius of a bound orbit. In region M3 
the potential has a local minimum and a local maximum and £ intersects V c g three times. The two intersection 
points at smaller p correspond to the minimal and maximal radius of a bound orbit, while the intersection point 
at larger p corresponds to the minimal radius of an escape orbit. Finally, in region M4 £ intersects V e g once. 
This intersection point corresponds to the minimal radius of an escape orbit. As already mentioned the form of 
the effective potential will change with 7 and j3. This leads then to a change of the p-v-pl&ne which is indicated 
in Fig|6ljb)-(e), where we demonstrate how the regions M1-M4 change when changing either 7 or /3. For fixed 7 
the regions M2 and M3 in which bound orbits exist get smaller when increasing /3 and in fact - as stated above 
- disappear completely for (3 > 2. This is seen in Fig|BJb)-(c). Moreover, for fixed (3 < 2 the regions M2 and 
M3 get bigger when increasing the gravitational coupling 7. This is shown in Fig[6jd)-(e). 

Note that our results for p z ^ are qualitatively similar, this is why we don't discuss them in detail here. 




Figure 4: The effective potential for the motion of a massless test particle (e = 0) in the space-time of an 
Abelian-Higgs string with f3 = 0.50 and 7 = 0.40 for different values of E. Here p z — 0.100 and L z = 0.015. 
The red horizontal line represents the value of £ . 
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Figure 5: The effective potential for the motion of a massive test particle (e = 1) in the space-time of an 
Abelian-Higgs string with j3 — 0.50 and 7 = 0.40 for different values of p z and L z . Here E = 1.060 and 
L z = 0.015. The black horizontal line represents the value of £ . 



The main difference between p z = and p z 7^ is that the particles move on 3-dimensional orbits in the latter 
case, while their motion is restricted to the x-y-plane in the former. 
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(a) /i-i'-plane for p = 0.5 and ■y = 0.4 




(b) p = 1.9, 7 = 0.4 



(c) = 1.5, 7 = 0.4 





(d) /3 = 0.5, 7 = 0.2 



(e) P = 0.5, 7 = 0.6 



Figure 6: (a): The different regions in the /i-^-plane for massive test particles (e — 1) with p z = 0. Bound 
orbits are possible in regions M2 and M3, while in Ml and M4 only escape orbits exist (see also FigJ7]for more 
details). The blue vertical line indicates the value of the energy at which the transition from region M3 to 
region M2 takes place, (b)-(e): The change of the regions M1-M4 when changing the gravitational coupling 7 
and the Higgs to gauge boson mass ratio /3, respectively. Note that for j3 > 2 no bound orbits exist. 
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Figure 7: The effective potential V c g(p) (blue) as well as (E 2 — e)/2 (red) are given for the different regions 
M1-M4 of the /x-z/-plot (see also FigJBJ). Note that for particle motion to be possible we need to require 
(E 2 — e)/2 > V^ff (p). Hence in M2 there is a bound orbit, in M3 a bound and an escape orbit, while in M4 
there is an escape orbit. 
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3.2 Examples of geodesies 

3.2.1 Massive particles e = 1 

In FiglSHTTl we show how bound and escape orbits of massive test particles (e = 1) change when changing 7 
and /?, respectively. In order to understand how the test particle moves we also indicate the radius of the scalar 
core (red) and of the magnetic flux tube core (blue). Since we measure the distance p in units of -M\y/v2 the 
scalar core radius is given by pn ~ W? an< ^ the magnetic flux tube core is given by /?w ~ l/v2- 

From Fig J8] it is evident that the maximal radius of a bound orbit decreases strongly with increasing 7. While 
it is much larger than pu for 7 = 0.36, it becomes comparable to pu for 7 = 0.42 and even smaller than pn 
for 7 = 0.48. This is related to the increased curvature of space-time for increasing 7. Moreover, the minimal 
radius of the escape orbit is always inside the magnetic flux tube core. Hence, the particle does not only move 
in the exterior (vacuum) region of the cosmic string, but can enter into the scalar and flux tube cores. 

The change of a bound orbit with (3 is shown in Fig{9j The maximal radius of the bound orbit increases with 
increasing (3, while the radius of the scalar core pn decreases at the same time. This leads to the observation 
that for small (3 the particle moves basically only inside the scalar core, while for increasing (3 it moves further 
and further away from the cosmic string. 

In Fig llOl we show how an escape orbit changes when changing the gravitational coupling 7. For 7 = 0.15 
and 7 = 0.3 the particle arrives from infinity, gets deflected by the cosmic string and moves again to infinity. 
The larger 7 the bigger is the change in the direction of the motion of the particle. This is not suprising since 
the deficit angle A increases with increasing 7. For 7 = 0.45 we observe a new phenomenon. The particle 
arrives from infinity, encircles the cosmic string and then moves again to infinity. This is new as compared to 
a space-time of an infinitely thin string. However, this type of encirclement has already been observed in the 
space-time of a Schwarzschild black hole and Kerr black hole, respectively pierced by an infinitely thin cosmic 
string for sufficiently large deficit angle [Til US] • Comparing these results with those given in Fig{TJJ we observe 
that the encirclement disappears when increasing (3. While it is still present for (3 = 0.5, it has disappeared for 
/3 = 2 and /3 = 5.3. 

3.2.2 Massless particles e = 

It seems that massless particles can only move on escape orbits, i.e. bound orbits as the ones shown for massive 
particles are not possible. We show the change of an escape orbit with the gravitational coupling 7 in FigfT2l 
For sufficiently large 7 we observe that the massless test particle encircles the cosmic string before moving again 
to infinity. This is completely new as compared to the space-time of an infinitely thin cosmic string. The change 
of the orbit with /3 is qualitatively similar to that of a massive particle. This is why we don't present it here. 

3.3 Observables 

3.3.1 Perihelion shift 

We can calculate the perihelion shift of a planar (p z = 0) bound orbit of a massive test particle (e = 1) by using 
([26]). The perihelion shift then reads 

Pmax 

8^ = 2 f ^ TT77-2-, (32) 

Pmin ^\P) \^ N( P ) 2 Hfi) 2 1 

where p m i n and p max are the minimal and maximal radius of the bound orbit. Our results for a particle with 
E = 1.03, L z = 0.1 and p z = are shown in Fig |13[ where the value of Sip is given in dependence on the 
gravitational coupling 7 for different values of (3. For this particular case, the perihelion shift is negative which 
means that the particle moves from p m - m to p ma , x and back to p m in under an angle of less than 27r. Though in 
most space times the perihelion shift is positive it is not too suprising that it can become negative in our case 
since the space-time is conical. Typically, the perihelion shift is negative for £ small, i.e. small values of the 
energy, while it becomes positive for larger values of £ . 
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(a) 7 = 0.36 
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(b) 7 = 0.42 
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(c) 7 = 0.48 

Figure 8: The change of the bound orbit of a massive test particle (e = 1) with E — 1.083, L\ — 0.025 and 
p z = 0.02 due to the change of the gravitational coupling 7. Here, we have chosen (3 — 0.5. The red and blue 
circle indicates the scalar field core and the magnetic flux tube core of the Abelian-Higgs string, respectively. 
The dotted lines indicate the minimal and maximal p of the bound orbit. In the 3-dimensional plots the red 
line indicates the string axis. 
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Figure 9: The change of the bound orbit of a massive test particle (e = 1) with E — 1.083, L\ — 0.025 and 
p z = 0.0205 due to the change of the ratio between Higgs and gauge boson mass /?. Here, we have chosen 
7 = 0.48. The red and blue circle indicates the scalar field core and the magnetic flux tube core of the Abelian- 
Higgs string, respectively. The dotted lines indicate the minimal and maximal p of the bound orbit. In the 
3-dimcnsional plots the red line indicates the string axis. 
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Figure 10: The change of the escape orbit of a massive test particle (e = 1) with E = 1.034, L\ = 0.02 and 
p z = 0.002 due to the change of the gravitational coupling 7. Here, we have chosen j3 = 1.25. The red and blue 
circle indicates the scalar held core and the magnetic flux tube core of the Abelian-Higgs string, respectively. 
In the 3-dimensional plots the red line indicates the string axis. 
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(b) /3 = 2.00 
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Figure 11: The change of the escape orbit of a massive test particle (e = 1) with E = 1.037, L\ — 0.015 and 
p z = 0.03 due to the change of the ratio between Higgs and gauge boson mass /3. Here, we have chosen 7 = 0.3. 
The red and blue circle indicates the scalar field core and the magnetic flux tube core of the Abelian-Higgs 
string, respectively. In the 3-dimensional plots the red line indicates the string axis. 
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Figure 12: The change of the escape orbit of a massless test particle (e = 0) with E — 0.075, L 2 Z — 0.03 and 
p z — 0.05 due to the change of the gravitational coupling 7. Here, we have chosen (3 — 0.5. The red and blue 
circle indicates the scalar held core and the magnetic flux tube core of the Abelian-Higgs string, respectively. 
In the 3-dimensional plots the red line indicates the string axis. 
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Figure 13: The dependence of the perihelion shift Stp on the gravitational coupling 7 for three different ratios 
/3 of the Higgs to gauge boson mass. Here, we have chosen e = 1, L z = 0.1, E = 1.03 and p z = 0. 

We observe that the value of the perihelion shift decreases for increasing 7. Moreover, increasing the ratio 
/3 between Higgs and gauge boson mass increases the value of the perihelion shift. 

3.3.2 Light deflection 



L =0.600, E = 1.400, P =0.000 

z z 




0.3 0.32 0.34 0.36 0.38 



Figure 14: The dependence of the light deflection 8(p on the gravitational coupling 7 for three different ratios j3 
of the Higgs to gauge boson mass. Here, we have chosen e = 0, L z = 0.6, E = 1.4 and p z = 0. 
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The deflection of light by a cosmic string can be calculated by using ([25]) for a planar (p z =0) escape orbit 
of a massless test particle (s = 0). The light deflection then reads 

oo 

r v = 2 I Lzdp , %1/a -ir, (33) 



^(P) 2 ( 






where p m i n is the minimal radius of the orbit. 

Our numerical results are shown in Fig |141 where we give the value of 5tp in dependence on the gravitational 
coupling 7 for three different values of /3. Apparently, light deflection becomes stronger when increasing 7 which 
results from the increase of the deficit angle. Moreover, increasing the ratio between Higgs and gauge boson 
mass increases the light deflection. 

4 Conclusions 

In this paper, we have studied the geodesic motion of massive and massless test particles in the space-time of an 
Abelian-Higgs string. We find that the existence of bound orbits, i.e. orbits on which a particle moves between 
a finite maximal and finite minimal radius depends crucially on the choice of the ratio between the symmetry 
breaking scale and the Planck mass and the choice of the ratio between the Higgs and gauge boson mass. In this 
paper, we have only considered the case of small symmetry breaking scale (in comparison to the Planck mass), 
i.e. we have only taken space-times with deficit angle smaller than 2n into account. Moreover, we have not 
studied Melvin-type space-times which we believe are physically not relevant from a astrophysical/cosmological 
point of view. 

We observe that bound orbits are only possible if the test particle is massive and if the Higgs boson mass is 
smaller than the gauge boson mass or in other words if the magnetic flux tube core lies inside the scalar core 
of the cosmic string. Increasing either the ratio between the symmetry breaking scale and the Planck mass or 
the ratio between the Higgs and gauge boson mass the particles move closer and closer around the string core. 
In fact, they move inside both the scalar and flux tube core and are not restricted to the movement inside the 
vacuum region outside both cores. For massless particles, which can only move on escape orbits, we observe a 
new phenomenon as compared to the space-time of an infinitely thin cosmic string: the particles can encircle 
the string before moving again to infinity. 

The perihelion shift of bound orbits of massive particles can be either negative or positive depending on 
the particle's energy and decreases with increasing ratio between symmetry breaking scale and Planck mass. 
Moreover, it increases with increasing Higgs to gauge boson mass ratio. The light deflection by a cosmic string 
increases with both the ratio between symmetry breaking scale and Planck mass and the ratio between Higgs 
and gauge boson mass. Since one of the possible detections of cosmic strings would be by light deflection 
(i.e. gravitational lensing) our results can be used to compare possible observational data to make predictions 
about the symmetry breaking scale at which the cosmic string formed as well as about the ratio between the 
corresponding Higgs and gauge boson mass of the underlying field theory. 

It would also be interesting to study geodesic motion in the space-time of a semilocal string [251 [UJ HH] , 
which is a solution of the electroweak model with gauge group SU(2)x U(l) in the limit where the Weinberg 
angle #vv = t/2. Moreover, field theoretical models describing so-called p-g-strings have been studied recently 
[29l [30] . p-g-strings are supersymmetric bound states of F- and D-strings and provide one of the examples 
of strings that might have formed in inflationary models resulting from string theory. Understanding how test 
particles move in the space-time of these strings would provide a further possibility to detect them. 
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